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Abstract 

o 

We study the entanglement creation in the low-energy scattering of two particles in three dimensions, for a general 

Ch . . . 

class of interaction potentials that are not required to be spherically symmetric. The incoming asymptotic state, before 

1—5 

the collision, is a product of two normalized Gaussian states. After the scattering the particles are entangled. We take 

,_i as a measure of the entanglement the purity of one of them. We provide a rigorous explicit computation, with error 

rS I bound, of the leading order of the purity at low-energy. The entanglement depends strongly in the difference of the 
I 

(— ( masses. It takes its minimum when the masses are equal, and it increases rapidly with the difference of the masses. 

a . . 

^ It is quite remarkable that the anisotropy of the potential gives no contribution to the leading order of the purity, on 
I— —I spite of the fact that entanglement is a second order effect. 

(N 
> 

in 

O 1 Introduction 
O 

In this paper we study how entanglement is created in a scattering process. This topic has intrinsic interest. As is 
well known, scattering is a basic dynamical process that is essential across all areas of physics. Entanglement is a 

^ central notion of modern quantum theory, in particular, it is the fundamental resource for quantum information theory 
and quantum computation. It is a measure for quantum correlations between subsystems. In the case of bipartite 

d systems in pure states, entanglement is a measure of how far away from being a product state a pure state of the 
bipartite system is. Product states are called disentangled. It is now well understood that entanglement in a pure 
bipartite quantum state is equivalent to the degree of mixedness of each subsystem. See, for example [1], [2], [3]. 
Moreover, the study of entanglement creation in scattering is interesting for a many other reasons. For example, for 
the implementation of quantum information processes in physical systems where scattering is central to the dynamics, 
like ultracold atoms and solid state devices. Moreover, the study of entanglement in the scattering of particles requires 
quantum information theory with continuous variables and mixed continuous-discrete variables. See [3] for a review of 
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this topic. As scattering interactions are fundamental at all scales, and as there is a large variety of scattering systems, 
it is possible that scattering will provide a new perspective to quantum information theory. Finally, entanglement 
creation is important to the theory of scattering itself, because it poses new problems that can shed some new light 
and new points of view in the study of scattering processes. 

Actually, from the conceptual point of view scattering is perhaps the simplest way to entangle two particles. Before 
the scattering, in the incoming state, the two particles are in a pure product state where they are uncorrelated. As 
they approach each other they become entangled by sharing quantum information between them. After the scattering, 
when they are far apart from each other, they remain entangled in the outgoing asymptotic state, that is not a product 
state anymore. 

We take as measure of entanglement of a pure state the purity of one of the particles, that is to say, the trace 
of the square of the reduced density matrix of one of the particles, that is obtained by taking the trace on the other 
particle of the density matrix of the pure state. The purity of a product state is one. 

We consider two spinless particles in three dimensions with the interaction given by a general potential that is not 
required to be spherically symmetric. Initially the particles are in an incoming asymptotic state that is a product of 
two Gaussian states. After the scattering the particles are in an outgoing asymptotic state that, as mentioned above, is 
not a product state and our problem is to determine the loss of purity of one of the particles, due to the entanglement 
with the other that is produced by the scattering process. 

The Hilbert space of states for the two particles in the configuration representation is 'H := L^(IR^). The Schrodinger 
equation is 

d 

«^^¥'(Xl,X2) = ff</7(xi.X2), (1.1) 

ot 



where the Hamiltonian is given by 



with Hq is the free Hamiltonian, 



H = Ho + V{^i-^2), (1.2) 



Ho:=- — A,-—A„ (1.3) 
2mi 27712 

where h is Planck's constant, mj,j = 1, 2, are, respectively, the mass of particle one and two, and Aj, the Laplacian in 
the coordinates Xj, j — 1, 2. The potential of interaction is multiplication by a real-valued function, V(x), defined for 
X G M'^. As usual, we assume that the interaction depends on the difference of the coordinates xi — X2, but no spherical 
symmetry is supposed. We assume that V satisfies mild assumptions on its regularity and its decay at infinity. See 



Assumption 2.1 in Section 2. For example, V{x) satisfies Assumption 2.1 if there are constants R,C > such that. 



'|x|<fl 



|\/(x)pdx < oo. 



(1.4) 



2 



and 



inx)i<c(i + |xi: 



-0' 



for |x| > R, 



(1.5) 



for some > /3, with (3 as in Assumption 2.1 Note that j3 controls the decay rate of the potential at infinity. Remark 



that (1.4) allows for Coulomb local singularities. We also suppose that at zero energy there is neither an eigenvalue 



nor a resonance (half-bound state), what generically is true. See Section 2. 

We work in the center-of-mass frame and we consider an incoming asymptotic state that is a product of two 
normalized Gaussian states, given in the momentum representation by, 



'/5in,po(Pl,P2) := (^po(Pl)'^-po(P2), 



(1.6) 



with. 



<Ppo(Pl) 



1 



= -(pi-Po)V2o-^ 



(a27r)3/4 

where p^, i = 1,2 are, respectively, the momentum of particles one and two. 



(1.7) 



In the state (1.6) particle one has mean momentum po and particle two has mean momentum — Pq. The variance 



of the momentum distribution of both particles is a. We assume that the scattering takes place at the origin at time 



zero, and for this reason the average position of both particles is zero in the incoming asymptotic state (1.6). After 
the scattering process is over the two particles are in the outgoing asymptotic state, (^out.po: given by 



'/fout,po(Pi,P2) := (5(pV2m)v3i„^pJ (pi,P2), (1.8) 

where p :— — ^ — Pi — — — Pt is the relative momentum, m := toi m-^/im^ + mo) is the reduced mass, and 
5(p2/2m) is the scattering matrix for the relative motion. 

The purity of (y5out,po is given by, 

^((Pout.po) = / rfPlc'pirfP2rfP2¥'out,po(Pl,P2)¥'out,po(p'l,P2)¥'out,po(pi,P2)¥'out,po(Pl,P2)- (1-9) 

Since the relative momentum, p, depends on pi and on p2, 'Pout,po is no longer a product state and it has purity 
smaller than one, what means that entanglement between the two particles has been created by the scattering process. 



Observe that in the state (1.6) the mean relative momentum of the particles is equal to po. 



Note that to be in the low-energy regime we need that the mean relative momentum po be small, but also that the 
variance a be small, because if a is large, the incoming asymptotic state ^in,-pa ^il^ have a big probability of having 
large momentum, even if the mean relative momentum po is small. 

We denote by (^in the incoming asymptotic state with mean relative momentum po = 0, and we designate, 
<Pout := 5(p2/2m)(/3i„. 
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We denote by 



TO,; 



TOi + m2 

the fraction of the mass of the i particle to the total mass. 



,i = 1,2, 



In Theorems |3.2| and |3.4| in Section 3 we give a rigorous proof of the following results on the leading order of the 
purity at low energy. 



^(V'out.po) = Viifiout) + O(|po/;i|), as |po/;i| ^ 0, 



(1.10) 



where O(|po/fi.|) is uniform on cr, for a in bounded sets. Furthermore, with /3 as in Asumption 2.1 (recall that 
controls the decay rate of the potential at infinity) , 



Viifiout) = 1 



(1.11) 



with Co the scattering length that is defined in (2.24 1 and where the entanglement coefficient £(/ii) is given by, 
16 4 



f (Ml) := 



ll + (2m-l) 



V 3/2 



- 1 



^ (1 + (2/.1 - 1)2) (2^1-1)2 ^1 + (2mi - 1)2 



8J{fii,l - fii) - 8J(1 - fii^fii) 



with 



Ji^J■l,^J'2) := ^ /rfq2 [ /dqi|A*2qi - Miq2| Exp[-1(m? + M|)(qi + q2)^ - (Ai2qi - Miq2)^ - q?/2] 

sinh[(/ii - /i2) |qi + q2| |M2qi - Miq2|] 

(mi - M2) |qi + q2| lM2qi - Miq2| 

In the appendix we explicitly evaluate J(l/2, 1/2) and J(l, 0), 



J(l/2,l/2) = ^ + - 



27 



3 arctan 



= 0.663497, 



(1.12) 



(1.13) 



(1.14) 



J(1,0) = 2(1 + ^ - ^2) = 0.32627. 
v3 

For /ii G [0, 1] \ {1/2, 1} we compute J(/ii, 1 — /ii) numerically using Gaussian quadratures. 



(1.15) 



Observe that £{ni) — £{1 — /ii), as it should be, because Viipout) is invariant under the exchange of particles one 
and two. 



Note that there is no term of order cr/h in (1.11 1. Actually, the terms of order a/h cancel each other because of 



the unitarity of the scattering matrix. This shows that for low energy the entanglement is a second order effect. 

The scattering length co is a measure of the strength of the interaction. As is well known, and can be seen in 
Theorem |2.2| in Section 2, at first order for low energy the scattering is isotropic and the total cross section, that is 
given by 47rco, is determined by the scattering length, cq. However, the effects of the anisotropy of the potential appear 
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at second order. It is quite remarkable that these effects give no contribution to the evaluation of the leading order of 



the purity. It follows from this that the leading order of the entanglement for low energy (1.11 1 is determined by the 
scattering length, cq, and that the anisotropy of the potential plays no role, on spite of the fact that entanglement is 
a second order effect, what is surprising. 

We see from Table 1 and Figure 1 that the entanglement coefficient depends strongly in the difference of the masses. 
It takes its minimun for /ii — 0.5, when the masses are equal, and it increases rapidly with the difference of the masses, 
as /ii tends to one. This shows that, if the scattering length is fixed, the entanglement takes its minimum when the 
masses are equal and that it strongly increases with the differences of the masses. This is indeed a remarkable result. 
Suppose that we consider different pairs of particles that interact in the same way at low energy, in the sense that, to 
leading order, they have the same total scattering cross section, i.e., such that the scattering length, cq, is the same 
for all the pairs. Moreover, suppose that the total mass, M, of the pairs is keep fixed, but that the individual masses, 
mi,rn2 of the particles are different in each pair. Our results show that, under these conditions, over four times more 
entanglement is produced by increasing the difference of the masses of the particles in the pairs. In practical terms, 
this means that in experimental devices to produce entanglement by scattering processes it is advantageous to use 
particles with a large mass difference. 

This fact can be understood in a physically intuitive way as follows: in the scattering of a light particle with a very 
heavy one, the trajectory of the light particle will be strongly changed, with a large exchange of quantum information 
between the particles, leading to a large entanglement creation. 

Note that in the scattering of a particle with a large mass and a particle with a small mass we can assume that 
the trajectory of the large particle is not affected by the interaction, i.e. that, to a good approximation, it follows a 
free trajectory, and that the small particle feels a (external) interaction potential centered in the position of the large 
particle. However, the trajectory of the small particle will be strongly affected by the interaction, what will produce 
exchange of information between the particles, leading to the creation of entanglement between them. To evaluate 
this entanglement it is, however, necessary to take into account the degrees of freedom of both particles, as we do to 
compute the purity. 

In the paper [4j a similar problem is considered in the case of equal masses and spherically symmetric potentials. 
They give an approximate expression for the leading order of the purity in the case of a Gaussian incoming wave 
packet that is very narrow in momentum space. 

The generation of entanglement in scattering processes has been previously considered in one dimension, mainly 
for potentials with explicit solution. See [5], [6], and the references quoted there. Moreover, [7], [8], [9], and the 
references quoted there, consider a system of heavy and light particles. They study the asymptotic dynamics and the 
decoherence produced on the heavy particles by the scattering with light particles in the limit of small mass ratio, 
what is different from our problem. The loss of quantum coherence induced on heavy particles by the interaction with 
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light ones has attracted much interest. See for example [TU], and [llj . 

The paper is organized as follows. In Section 2 we define the wave and scattering operators, the scattering matrix 
and we consider its low-energy behavior. In Section 3 we prove our results in the creation of entanglement. In Section 
4 we give our conclusions. In the Appendix we compute integrals that we need in Section 3. Along the paper we 
denote by C a generic positive constant that does not necessarily have the same value in different appearances. 

2 Low-Energy Scattering 

We consider the scattering of two spinless particles in three dimensions. We find it convenient to use the time-dependent 
formalism of scattering theory. See, for example, UM [T51 [TBI 117) . 



The Hilbert space of states in the configuration representation is "H ;= L (IR°). The Schrodinger equation is 



where the Hamiltonian is given by 



d 

ih—ip{xi,X2) = Hip{xi.X2), (2.1) 



H = Ho + V{xi~X2)- (2.2) 



The operator Ha is the free Hamiltonian, 

H,:=- — A,-—A2, (2.3) 
2mi 2ni2 

with h Planck's constant, mj,j = 1,2, respectively, the mass of particle one and two, and Aj the Laplacian in the 
coordinates Xj,j = 1,2. The potential of interaction is multiplication by a real-valued function, V(x), defined for 
x e M'^. As usual, we assume that the interaction depends on the difference of the coordinates Xi — X2, but no spherical 
symmetry is supposed. V satisfies the following condition. 

ASSUMPTION 2.1. For some (3 > 0, {1 + \x\)^V{x) is a compact operator from the Sobolev space into the 
Sobolev space . 

Below we will assume that /3 > 5 or that /3 > 7. For the definition of Sobolev's spaces see [H]. Conditions for 



Assumption 2.1 to hold are well know [TB], [12]. For example, if (1.4, 1.5) are satisfied. 

Under this condition H is defined as the quadratic form sum of Hq and V and it is a self-adjoint operator. 
The wave operators are defined as 

As is well known, under our condition the wave operators exist and are asymptotically complete, i.e. their ranges 
coincide with the absolutely continuous subspace of H. Moreover, the scattering operator, 

S^W^W., (2.4) 



is unitary. 

Before the scattering, when the particles are far apart from each other and the interaction is weak, the dynamics 
of the system is well approximated by an incoming solution to the free Schrodinger equation with the potential set to 
zero, 



where the incoming asymptotic state (p_ is the Cauchy data at time zero of the incoming solution to the free Schrodinger 
equation. When the particles are close to each other, and the potential is strong, the dynamics of the system is given 
by the solution to the Schrodinger equation, 

e-'i"W-ip-, (2.5) 
that is asymptotic to the incoming solution to the free Schrodinger equation as t — >■ — oo, 

e -'i"" . 



lim 



0. 



After the scattering, for large positive times, the particles again are far away from each other and the dynamics of the 
system is well approximated by the outgoing solution to the free Schrodinger equation 



that is asymptotic to the solution to the Schrodinger equation (2.5) as t — >■ oo 



lim 



The outgoing asymptotic state is the Cauchy data at time zero of the outgoing solution to the free Schrodinger 
equation; ip^ := W-cp-. It is given by the scattering operator, ip^ — Sip-. 



= 0. 



As usual, we consider the center-of-mass and relative distance coordinates, 

miXl + 7712X2 

Xcm ■ J 

nil + 1712 
X := xi - X2. 



(2.6) 
(2.7) 



The state space H factorizes under this change of coordinates as, 



H — Hem (8) Hid, 



(2.8) 



Where Hem — L'^{^^),Hic\ '■— L^(IR'^) are, respectively, the state spaces for the center-of-mass motion and the relative 
motion. Since the interaction depends only on x, the Hamiltonian and the wave and scattering operators decompose 
under the tensor product structure and, in particular, we have that 



(2.9) 



where /cm i« the identity operator on Hem and S'rci is the scattering operator for the relative motion, that is defined 
as follows. The Hamiltonian for the relative motion is given by, 

ifrei:=-|^A^ + y(x), (2.10) 
2m 

where m is the reduced mass, 

mi 1712 

m:= , (2.11) 

TOi + m2 

and Ax is the Laplacian in the x coordinate. The free relative Hamiltonian is, 

Ho,rel ■■= -|^^x. (2.12) 

The relative wave operators are defined as, 

W±,rel := S-lim±^ e'S^-' g-iiHo,re/_ (2.13) 

The relative scattering operator, 

5,el = VF;,,elW^-,rel, (2.14) 

is a unitary operator on Hrei- 



We denote by H := X^(M^) the state space in the momentum representation. The momentum of the particles one 
and two are, respectively, pi,P2. We define the Fourier transform as an unitary operator from H onto H, 



^¥'(Pi,P2):-7;^ / e-H(P-W--=)<p(xi,X2). 



(2.15) 



It is also convenient to take as coordinates in the momentum representation the momentum of the center of mass and 
the relative momentum, 

Pcm:=Pl+P2, (2.16) 

p^^m2m^miP2_ 
mi + m2 

The state space in the momentum representation also factorizes as a tensor product, 

U = Ucu.®UrA, (2.18) 

where i-Lcm = -t'^(R^),^rei := i'^(M') are, respectively, the state spaces in the momentum representation for the 
center-of-mass motion and the relative motion. 

The scattering operator in the momentum representation, 

S:=J'SJ^-\ (2.19) 

decomposes as, 

5 = /em 5,el, (2.20) 



where S^oi is the scattering operator for the relative motion in the momentum representation, 

S'rcl ■— J'rcl Srcl J'^el ^ (2-21) 

where J-'rci is the Fourier transform in the relative coordinate, 

J-rei^(p) — ^ / e-*P-Xx). (2.22) 

We denote by the unit sphere in M.^. 

As 5rci commutes with -ffo,rci (energy conservation) we have 

(5,.ei¥>) (p) = (5(pV2m)(p) (p), (2.23) 

where the scattering matrix, S{E), is a unitary operator in i^(§^) for each E e (0,oo). Note that the scattering 
matrix defined in the time-dependent framework coincides with the scattering matrix defined in the stationary theory 
by means of the solutions to the Lippmann-Schwinger equations. 

The following theorem has been proved by Kato and Jensen [TO]. Note that they consider the case h = 1, to = 1/2 
but the general case is easily obtained by an elementary argument. A zero energy resonance (half-bound state) is a 
solution to -ffre;</5 = that decays at infinity, but that is not in L'^{M.^). See [TO] for a precise definition. For generic 
potentials V there is neither a resonance nor an eigenvalue at zero for H^^i . That is to say, if we consider the potential 
AV^ with a coupling constant A, zero can be a resonance and/or an eigenvalue for at most a finite or denumerable set 
of A's without any finite accumulation point. 

The scattering length is defined as, 

1 /2to , 2to , 1 \ , , 

co:=- [^-j_^V{l + Go^Vr\iy (2.24) 

where (., .) is the scalar product in M.^, 1 designates the fimction identically equal to one and Go is the operator 
with integral kernel the Green's function at zero energy, 

Go(x,y):=— ^ |,x,yeM3^ (2.25) 

47r|x-y| 

The operator (1 + Gq^V^) is invertible because zero is neither an eigenvalue nor a resonance for H^^i- We define the 
scattering length, co, with the opposite sign to the one used in ,19J . in order that it coincides with the definition used 
in the physics literature [TOIIIT]. Furthermore, 



ro(j^) := ^,J^e§', (2.26) 



1 I 2m,, / „ 2to. ^ ^ 



I 7^ (l + Go^V^) l,x.H-e§^. (2.27) 



We denote by B (i^(§^)) the Banach space of all bounded linear operators on L^(§^). 



THEOREM 2.2. (Kato and Jensen 119^ ) Suposse that Assumption 2.1 is satisfied and that at zero H^^i has neither 
a resonance (half-bound state) nor an eigenvalue. Then, If jS > 5, in the norm o/S(L^(§^)) we have for \p/h\ — >■ 
the expansion, 

5(pV2m) = J + z\p/h\ - \p/hf S° + o{\p/hf), (2.28) 
where I is the identity operator on L^(§^), 

E; :--2co(-,ro) 1^, (2.29) 

and 

:= 2cl (., Yo) Yo + (•, Fi) Fq - (•, Ya) Y,. (2.30) 
Furthermore, if P > 7, o(|p//ip) can he replaced by 0{\p/h\^). 

Note that Yi = if is spherically symmetric. We see that, as is well known, the leading order at low energy of 
iS(p^/2m) — / is given by the scattering length, i.e. in leading order the scattering is isotropic. The anisotropic effects 
appear at second order. 

3 Entanglement Creation 

Consider a pure state of the two-particle system given in the momentum representation by the wave function ip{pi, P2)- 
Let us denote by p{ip) the one-particle reduced density matrix with integral kernel, 



p{'P){Pi,p'i) </5(Pi,P2)¥'(pi,P2) dp2, 



and by 'P{(p) the purity, 

'P((^) := Tr(p2) = / dpidpidp2dp^(/3(pi,P2)v3(pi,P2)</5(pi,P2)'y5(Pi,P2)- (3-1) 

The purity is an entanglement measure that is closely related to the Renyi entropy of order 2,— In Tr(/3^) [3l [20l |2T] . 
It is trivially related to the linear entropy, Sl, SisSL — i — V. It satisfies < < 1 if <y9 is normalized to one. 
Furthermore, it is equal to one for a product state, ip = ipi{pi) (y92(p2). The purity is an entanglement measure that 
is convenient for the study of entanglement creation in scattering processes because it can be directly computed in 
terms of the scattering matrix. 

We work in the center-of-mass frame and we consider an incoming asymptotic state that is a product of two 
normalized Gaussian wave functions, 

'Pin,po(Pl,P2) := (^po(Pl)¥'-po(P2), (3.2) 



where 



V'poCpO-P^-'^^^-^-^^/'^^ (3-3) 
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In the incoming asymptotic state (3.2) particle one has mean momentum pg and particle two has mean momentum 
— Pq. The variance of the momentum distribution of both particles is a. We assume that the scattering takes place at 
the origin at time zero, and for this reason the average position of both particles is zero in the incoming asymptotic 



state (3.2). Note that by (2.17) the mean value of the relative momentum in the state (3.2) is equal to pp. 
Since ^in.pa is a product state its purity is one, 

^(^in.po) = 1- (3.4) 

After the scattering process is over the two particles are in the outgoing asymptotic state, </'out,poi given by 

<y5out,po(Pl,P2) (5(pV2TO)¥'i„,po) (Pl,P2)- (3.5) 

Since the relative momentum, p, depends on pi and on p2, •fout.-pa is no longer a product state and it has purity 
smaller than one, what means that entanglement between the two particles has been created by the scattering process. 

We will rigorously compute the leading order of the purity of lySout.po "in a quantitative way- in the low-energy limit 
for the relative motion. Note that to be in the low-energy regime we need that the mean relative momentum po be 
small, but also that the variance a be small, because if a is large the incoming asymptotic state </'in,po "^iU have a big 
probability of having large momentum, even if the mean relative momentum po is small. 

We first introduce some notations that we need. 

We denote by (pin the incoming asymptotic state with mean value of the relative momentum zero, 

<^in(Pi, P2) := ¥'(Pi) </'(P2), (3.6) 

where, 

^(P):= (^^4)^e-PV2.^,peM^ (3.7) 
and by (pout the outgoing asymptotic state with incoming asymptotic state (/7in, 

<<3out(Pi,P2) (5(pV2m)<^i„) (pi,P2). (3.8) 

We define. 



V'qo(q):= (^e-^'i-'i'') /^qeE^ (3.9) 

V'(q):=^^e-'i'/2,qeM3, (3.10) 

^in,qo(qi,q2) ^Aqo(qi) V-qo(q2), (3-11) 

V'i„(qi,q2) := V^(qi)^(q2). (3.12) 



We prepare the following proposition that we need later. 
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PROPOSITION 3.1. 



Proof: We denote qo := po/a. Then, 



llv'in.po -^inW < C'min{|po|/f7, 1}, 

llp(¥'m,po - <y'm)ll < C\po\. 
ll'Pin.po - t^inll = llV'in.qo " "^inH ■ 



Assume first that |qo| < 1. We have that, 



Moreover, 



g-qo+(qi-q2)-qo _ i 



-qo+(qi-q2)-qo 



ds 



< 



Jqo|^ + |qo|(|q.| + |q.|)(|q^|2^|q^|(|q^| + |q^|)). 



(3.13) 
(3.14) 

(3.15) 

(3.16) 

(3.17) 
(3.18) 



It follows from ( |3.15|3.16[ and |3.18[ ) that ( |3.13| holds for qo| < 1. In the case |qo| > 1 the estimate is immediate 
because, 

||</'in,po - <Pin|| < 2. 



Note that by (2.17), 



IpI < |pi| + |p2| 



Then, if |qo| < 1 as in the proof of (3.13) we prove that, 

||p(¥'in,po - Vin)\\ < C\po\. 

If IqoI > 1 we estimate as follows, 

||p(<y5in,po - -^in)!! < |||Pl| (Vin.po " -^in)!! + |||P2| (<y5i„,po - 'fiiiM 



(3.19) 



(3.20) 



Furthermore, 



|Pl| (tPin^po - Vin)\\ < IIIPI - Pol Vin,po|| + |||Po| tPin,poll + ll|Pl| finW < 

c^lllqi - qol V'in,qoll + IpoI + O-|||qi|'0in|| < C|po|. 



(3.21) 



In the last inequality we used that a < |po|. In the same way we prove that. 



P2I (V'in.po - Vin)\\ < C\po\. 



(3.22) 



By (3.20 3.21 3.22) we have that. 



lp(<^in,po - Vin)\\ < C\po\. 



(3.23) 
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Equation (13. 141) follows from (|3.19|) and (13.23 1 



Let us denote, 



□ 



(3.24) 



where / designates the identity operator on L (§ ). It follows from (2.281 and since ||5(p /2m)||e(i2(s2')') = 1, that 



|r(pV2m)| 



< c 



lp/^l 



B(L2(s^)) -^i + |p/;i|- 



Hence, 



|r(pV2m)(piJ| < C- llq^i, 



(3.25) 



(3.26) 



We designate, 



C{(pl,<l>2,(f>3,<f>i) / C^PlrfpirfP2rfP2'/'l(Pl.P2)02(pi,P2)'/'3(pi,P2)'/'4(Pl,P2)- (3.27) 



Note that. 



7'(</)) = £(</), <^,0,(/.). 

It follows from the Schwarz inequality that, 

|>C(0i,02,03,04)| <n|^l||0j| 

The following theorem is our first low-energy estimate of the purity. 



(3.28) 



THEOREM 3.2. Suposse that Assumption 2.1 is satisfied and that at zero H^-ei has neither a resonance (half-bound 
state) nor an eigenvalue. Then, 



^(V'out.po) = Viipout) + O(|po Al), as |po/;i| ^ 0, 
where O{\po/h\) is uniform on a, for a in bounded sets. 



(3.29) 



Proof: Writing v?out,po as, 



<Pout,po := S{p^/2m)ipin,po = <y5in,po + T(p^/2m)(/3in,po 



and using (3.4), we see that we can write 7'(<Pout,po) as follows, 

7'(v?out,po) = 1 +7^(po), 

where 72.(po) is given by, 

A 

7?-(Po) := ^A(P0, V'l, ^^2, V'3,V'4), 

1=1 

for some integer A, and where each of the £i(po, "(Aii "02, "03 7 "04) is equal to, 

^,:(PO,'01,V'2,V'3, 04) = -C(-0i,02,03,V'4), 



(3.30) 



(3.31) 



(3.32) 
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where for some 1 < k < A, k oi the tpj are equal to 7~(p^/2TO)(^in and the remaining 4 — fc are equal to V'in.po- 



Similarly, 



with 



^Kut) = i + n{o), 



7^(0) :=^A(0,V'i,V'2>3>4) 

i=l 



Below we prove that, 



7^(po) = 7^(o) + o(|po/;i|), as |po/;i| ^ o, 

what proves the theorem in view of p.30|3.33D . 



(3.33) 



(3.34) 



(3.35) 



We proceed to prove (3.35). Without losing generality we can assume that. 



£l(pO, -01, "02, V'3, ^^4) = /:(T(p^/2m)(pin,po, </5i„,po, (^in,po, (^in^po)- 



(3.36) 



We have that. 



'Ci(po,V'i,'02,'!/'3,V'4) = C{T{p'^/2m)ip ) +£(r(pV2m)(</?i„,p„ - (^in),i^in.po,'^in,po,<^in,po)- 



(3.37) 



By (3.14 3.25 3.28 3.37), 



£l(pO,V'l,V'2,'03,'04) ^ ^T{p'^ /2m)ipin,ipin,po,'Pin,po,'Pin,po) + 0{\po/h\),a.s\po/fl\ 0. (3.38) 



In the same way, using (3.13 3.26|3.38 ), we prove that, 

^l(PO,'01,V'2,'03,'04) = .C(r(pV2TO)Vin,</3in,¥'in,po,<^in,po) + 0( |po//l| ) , aS |po//l| 0. 

Repeating this argument two more times we obtain that, 

£l(pO,-01,V'2,l/'3, V'4) = £(T(p^/2m)(/?in,(/Jin,<i5in,</?in,) + 0{\po/h\) = 

>Ci(O,V'i,'02,'03,V'4) +O(|po//i|),as|po/;i| -> 0. 

We prove in the same way that, 

£j-(po, V'i>2, V'3,'04) = £j(0, V'l, V'2,V'3,V'4) + O(|po/^|),2 <j< A,a.s\po/H\ 0. 



(3.39) 



Equation ( |3.35[ ) follows from ( |3.30|3.31|3.33|3.34[ [3T0|3.41 1 



(3.40) 



(3.41) 



□ 



We now compute the leading order of the purity of (pout- 
We denote, 



ri(pV2m) 5(pV2m) - / - i\p/h\ + \p/h\'' E°. 



(3.42) 
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It follows from Theorem 12.21 that 



|ri(pVm)| 



< 



\p/h\^o{l), if/? > 5, 



lB(i^(S^)) ^ 1 |p/;i|20(|p/fi|), if^ > 7, 

where o(l) and 0(|p//i|) are bounded functions of |p/^|, lim|p/ft|^o o(l) = and 0(|p//i|) < C|p//i| for |p/S| < 1 



(3.43) 



THEOREM 3.3. Suposse that Assumption 2.1 is satisfied and that at zero H-^ei has neither a resonance (half-bound 



state) nor an eigenvalue. Then, as a/h goes to zero, 



P(^out) + i\p/h\ - \p/h\^ s^] (pi„) 



o{\a/n\^), if/3>5, 
0{\a/h\^), if/3>7. 



(3.44) 



Proof: We write ip out as follows, 



'Pout = 'Pout,! + 7'i(pV2m)(^i, 



where, 

<^out,i := [/ + t\p/n\ S? - \p/h\^ E^] cpi, 
Using this decomposition we write V{(pout) as follows, 



(3.45) 



V{(pout) ^V{(p out.l 



(3.46) 



where 7?,((t) is given by. 



for some integer B, and where each of the Ci{a, ipi,ip2, "03) V'4) is equal to, 

where for some 1 < fc < 4, fc of the ipj are equal to (^out.i and the remaining 4 — fc are equal to 7i(p^/2m)(/3i, 



(3.47) 



(3.48) 



We proceed to prove that as cr/h goes to zero. 



Tl{a) = 



_( o{\a/h\^), if/3>5. 



0{\a/h\^), if/3>7. 



what proves the theorem in view of ( 3.461. 



Without any loss of generality we can assume that, 



(3.49) 



£b(0-, -01, "02, Tpa, 04) = >C((pout,l, </5out,li '/'out,!, 71(p^/2m)(^in). 



(3.50) 



By (3.28 3.43) We have that, 



r ( I I I I \ f o{\a/h\^), if/3>5, 
£5(^,01,02,03,04) = I o\\a/h\% if/3 > 7. 



We complete the proof of (3.49) estimating the remaining terms in (3.47) in the same way. 



(3.51) 
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We denote by 



Mi 



mi + m2 

the ratio of the mass of the i particle to the total mass. 



,1 = 1,2, 



It follows from ( |2.16|2.17p that 



and that, 



(cr27r)3/2 

By a simple computation using (3.45 3.53|3.55 ) we prove that. 



Pi = A^lPcm + P, 



P2 = M2Pcm - P, 



1^1+1^2 p2 P -^2)Pciii-P 



( [/ + t\p/h\ S? - |p/fi|2 (pi„) = 1 ^ (a/fi)2 (Pi(Vin) + P2(V^i„)) + asa/H ^ 0, 



where. 



Pl(Vin) = S|=lPlj(V'in), 



with 



'Pi,i(V'in) = '^qidq2dq3 |Ai2qi - A*iq2| lM2q3 - Aiiq2| (s?V'(qi,q2)) (s'i'-(/'(q3,q2))'0(qi,q3), 
^i,2(?/'in) = -2 y dqidq2dq3 |/i2qi - Miq2| lAi2qi - A'lqal (s^V'lqi, q2)) (s??/'(qi,q3)) '/'(q2,q3), 



^1,3 (V-h 



y" dqidqa |Ai2qi - /iiq2| (S?-0(qi, q2)) '0(qi,q2) 



and 



P2(V'in) = 4 y dqidqz |M2qi - Miq2p (SaV'lqi, q2)) V'(qi,q2)- 

Explicitly evaluating the integrals in (3.58), 3.59|3.60 1 using (3.55) and /i2 = 1 — Mi: prove that. 



7'i,i(V'i„) = -8^ J(/ii,l-Mi), 



7'i,2(^i„) = -8^ J(l-Mi,Mi), 
Pl,3(V'in) = 8^ (L(Mi,l-Mi)f , 



where. 



./(mi,M2) := ^ /c^q2 [ /rfqi|M2qi - Miq2| Exp[-i(^? + ^|)(qi + q2)2 - (^2qi - Miq2)^ - q?/2] 

sinh[(/xi - ^2) |qi + q2| lM2qi - Miq2|] 
(mi - M2) |qi + q2| lM2qi - Miq2| 
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and 



L(/ii,/^2) ;p3 /dqidq2|M2qi - ^^l<i2\ Exp[-(^f + ^|)(qi + qa)^ - 2(^2qi - A*iq2)^] 

sinh[(/ii - ^2) |qi + q2| lM2qi - Miq2|] 



Exp[-(Aii - Ai2)(qi + q2) • (/^2qi - Miq2)]- 



(mi - M2) |qi + q2| |M2qi - Miq2| 



Furthermore, 

where, 

^(^1,^2) 



2 2 

7'2(V'm)=8^7V(/ii,l-A'l), 



dqcm dq q^ Exp[-(^^ + fiDul^ - 2q^ - - ^l2)c^cm ■ q] 



smh[(^i - ^2) |qcm| |q[ 

(mi - A*2) |qcml|q| 



(3.66) 



(3.67) 



(3.68) 



Note that the second and the thhd term m the right-hand side of (2.30) give no contribution to A^(/ii,/i2) because as 
Yi(;/) is an odd function the integrals of these terms are zero. 



By ( |3l4l [3l6l [slTj [3!62pJ4l [3^ 



Vi^ont) - 1 - 8(coa/;i)2 {{L{ni, 1 - /xi))2 + iV(^i, 1 - /ii) - J(/ii, 1 - A^i) - J(l - + 

o(|a/;i|2), if/3>5, 
0(|a/;i|3), if/?>7. ^^-^^^ 



In the appendix we prove by expUcit computation that, 



L{^,^,l^^,,)^^ - (1 + (2/ii - 1)2) 



2N-1/2 



2(2a*i - 1)2 ^1 + (2^1 -1)2 
7V(l/2,l/2) = 3/4. 



(1 + (2^1 -1)2)^/^-1 



We denote by f (/ii) the entanglemement coefficient, 



£i^ll) := 8 [(L(/^i, 1 - Aii))2 + iV(^i, 1 - ^1) - J(mi, 1 - A^i) - ■/(! " • 



(3.70) 

(3.71) 
(3.72) 

(3.73) 



By ( |3.70|3.71D , 

£(/ii' - 
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(1 + (2mi- 1)2)^/^-1 



+ 



TT (1 + (2/ii - 1)2) (2^1-1)2 ^1 + (2^1 -1)2 

Thus, we have proven the following theorem. 



8J(Mi,l-Aii)- 87(1-^1,^^1)- (3.74) 



THEOREM 3.4. Suposse that Assumption 2.1 is satisfied and that at zero H-^ei has neither a resonance (half-bound 
state) nor an eigenvalue. Then, as a/h goes to zero, 



(3.75) 



where the entanglement coefficient £ (fii) is given by {3.74 ) 
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Proof: The theorem fohows from (3.69 3.731 



In the appendix we exphcitly evaluate J(l/2, 1/2), 



3 arctan 



n/3 



= 0.663497. 



By (3.74 3.761 for fii — 1/2, when the masses are equal, the entanglement coefficient is given by 

f (1/2) = 0.4770. 

We also explicitly evaluate in the appendix J(l, 0), 

J(1,0) = 2{l + ^-V2) = 0.32627. 
v3 

For /ii e [0, 1] \ {1/2, 1} we compute J(/ii, 1 — /ii) numerically using Gaussian quadratures. 
In Table 1 and in Figure 1 we give values of £{fJ-i) for 0.5 < /ii :— mi/{mi + TO2) < 1. 



□ 



(3.76) 



(3.77) 



(3.78) 



REMARK 3.5. Note that there is no term of order a/h in (3.75). Actually, the terms of order a/h cancel each 



other because of the complex conjugates in the definition of the purity in (3.1 1 and of the factor i in the second term 



in the right-hand side of (2.281 that is there because of the unitarity of the scattering matrix. This shows that for low 



energy the entanglement is a second order effect. 



REMARK 3.6. Remark that £{ni) = £{1 — /^i), what implies that the leading order in (3.75) is invariant under the 
change /^i f-s- 1 — /ii , as it should be, because 'P{(pout) is invariant under the exchange of particles one and two. 



REMARK 3.7. As we mentioned in Remark 3.5, at low energy the entanglement is a second order effect. As can 



be seen in Theorem 2.2 for low energy the scattering is isotropic at first order and it is determined by the scattering 
length, cq. However, the effects of the anisotropy of the potential appear at second order. It is quite remarkable 
that these effects give no contribution to the evaluation of N{^i, ^2), as mentioned above. It follows from this that 



the leading order of the entanglement for low energy (3.75) is determined by the scattering length, cq, and that the 



anisotropy of the potential plays no role, on spite of the fact that entanglement is a second order effect, what is 
surprising. 



4 Conclusions 



We considered the entanglement creation in the low-energy scattering of two particles of mass mi,rn2: in three 
dimensions with the interaction given by potentials that are not required to be spherically symmetric. Initially the 
particles are in a pure state that is a product of two normalized Gaussian states with the same variance, tr, and 
opposite mean momentum. The entanglement creation by the collision was measured by the purity, V ^ of one of the 
particles in the state after the collision. Before the collision the purity is one. 
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We gave a rigorous computation, with error bound, of the leading order of the purity, V , at low-energy. Namely, 
we proved that the leading order of the purity is given by 1 — (cocr//i)^£, where Cq is the scattering length and the 
entanglement coefficient £ depends only on the masses of the particles. 

We proved that the entanglement takes its minimum when the masses are equal and that it strongly increases with 
the differences of the masses. There is no term of order a /h in the leading order of the purity, what shows that for low 
energy the entanglement is a second order effect. As is well known, for low energy the effects of the anisotropy of the 
potential appear at second order. It was found that these effects give no contribution to the evaluation of the leading 
order of the purity and that the anisotropy of the potential plays no role, on spite of the fact that entanglement is a 
second order effect, what is surprising. 

5 Appendix 

For the reader's convenience we compute on this appendix the integrals that we need in Section II. 
We first state some elementary integrals that we need. 

e-'^^'dx^ a>0, (5.1) 

2 V a 



1 /tt 
4^V a' 



e dx ^ —\l a > 0, 



-ax —2hx 



dx^ ^|-e^^''', a>0, 
a 



-ax^-2hx rL gfcVa „ erf(&/\/^)), « > 0, 

2 V a 



where erf(x) is the error function. 



erf(a;) 



e ^ dy. 



(5.2) 
(5.3) 
(5.4) 

(5.5) 



Equation (5.2 1 follows integrating by parts using —(l/2a)^e = e a; and (5.1). Equations (5.3 5.4) follow 



from (5.1 ) changing the variable of integration io y = x -\- b/a. 



I (l + a;2)v 



(l + a;2) V(2 + a;2) 



dx = arctan 



arctan 



2 - V2 + a;2 



C. 



(5.6) 



We have that, 



^(erf(y^))^-- 



4 



dye 



Hence, 



1 - (crf(^/J))2 = 



4 /-i £--=^^"+1) 



|-(erf(^/^))2 = 

az TT Jq + 1 



dy. 
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It follows that, 



(erf(z))2 = 1 - - 



dy. 



(5.7) 



We first compute J(l/2, 1/2). By 



J(l/2,l/2) = 



47^9/2 



dq2e-'i^g(|q2|)' 



(5.8) 



where, 



5(|q2|) ■■= j (iqi|qi-q2|e dqi. 



(5.9) 



To evaluate (5.9 ) we take a system of coordinates where q2 = (|q2|, 0, 0), we do the change of coordinates (gi.i, 91.2, 91.3) 
" IQ2I7 '71,2, 91,3) and we compute the integral in spherical coordinates, to obtain. 



— |q2|^ poo 

5(|q2|) = / e-^^ p2 (^2|q.|p _ ^-2|q.|p)^ 

1^2 1 Jo 



After repeated integrations by parts using — (l/2a) J^e "-^ = e p and (5.4) we prove that, 



5(|q2|) 



2|q2| 



erf(|q2|)(l + 2|q2n+^e-l'i^l' 



Introducing (5.11) into (5.8) and passing to spherical coordinates we obtain, 



J(l/2,l/2) = — / -2.-"' 



Jo 



P e 



3/2 \2 

erf(p)(l + 2^2)4.776-" dp. 



2p 



After expanding the square in the right-hand side of (5.12 ), several integration by parts using — (l/2a) 



(5.10) 



(5.11) 



(5.12) 



-ap _ p — ap 



and (5.1 5.2 5.6 5.7 1 we obtain that 



J{l/2, 1/2) = \ + -{ ^ - Sarctan ( ^ ) ) = 0.663497. 



2 TT \ 4 



2-^3 



(5.13) 



We now compute J(1,0). By (3.65) 



•^(i'0) = 4^ / ^q2q2e-'''^ (Mq2))^ 



(5.14) 



where. 



, , ,2, 2, piqi+qal Iqal _ p-|qi+q2| |q2| 
h{ci2) = / dqie-^'i^+'i^) /'e-i?/2 , 

qi + q2 q2 



Changing the integration coordinate in (5.151 to Q = qi + q2 we obtain 



olQI|q2l _p-|QI|q2l 



=Q q2 



IQI|q2| 

Using spherical coordinates and doing the integration in the angular variables we get 

2ti 



Mq2) = ^e-^i^/^ / dpe-f^' 

|q2r Jo 



g2p|q2| ^ g-2p|q2| _ 2 



o-q2/2 



q2 



dpe-P 



o2p|q2| _ 



By ( [5l|531 ) 



Mq2) = ^e-'i^/' ie^^^i). 
|q2r 



(5.15) 



(5.16) 



(5.17) 



(5.18) 
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Introducing (5.18) into (5.14| and performing the remaining integrals with the aid of (5.1 ) we prove that 



J(1,0) = 2(1 + ^ - ^2) = 0.32627. 
v3 



(5.19) 



We proceed to compute i(/i, 1 — /ii). Using spherical coordinates and performing the integrals in the angular variables 
we prove that for fii ^ 112, 

4 



dX A e 



(5.20) 



By (|5.1 5.3l), and integrating by parts using —(l/2a)g^e — e ""P p we prove that, 



L(/ii,l-^i) = (1 + (2^1 -1)2)- 



-1/2 



(5.21) 



Finally, we compute iV(/ii, 1 — /ii). Using spherical coordinates in (3.68) and evaluating the integrals in the angular 
coordinates we obtain for /ii 7^ 112 that, 

4 



By ( |5.2|5.3[ ) 



7r(/ii - Ai2)2 7o 
iV(Mi,l-Mi) ^ 



Taking the limit as ^1 — > 1/2 we get. 
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Table 1: The Entanglement Coefficient £{ni) 



Hi := mil [mi 

0.5 

0.525 

0.55 

0.575 

0.6 

0.625 

0.65 

0.675 

0.7 

0.725 

0.75 

0.775 

0.8 

0.825 

0.85 

0.875 

0.9 

0.925 

0.95 

0.975 

1 



TO2) £{ni) 
0.4770 
0.4813 
0.4937 
0.5144 
0.5434 
0.5816 
0.6296 
0.6880 
0.7550 
0.8320 
0.9179 
1.0120 
1.1130 
1.2208 
1.3228 
1.4488 
1.5659 
1.6832 
1.8010 
1.9168 
2.0287 




